Renormalized mean-field theory of the neutron scattering in cuprate superconductors 
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The magnetic excitation spectrum of the t-t' -J-model is studied in mean-field theory and compared 
to inelastic neutron-scattering (INS) experiments on YBa2Cu306+ M (YBCO) and Bi2Sr2CaCu2C>8+<s 
(BSCCO) superconductors. Within the slave-particle formulation the dynamical spin response is 
calculated from a renormalized Fermi liquid with an effective interaction ~ J in the magnetic 
particle-hole channel. We obtain the so-called "41 meV resonance" at wave vector (it, n) as a 
collective spin-1 excitation in the d-wave superconducting state. It appears sharp (undamped), if 
the underlying Fermi surface is hole-like with a sufficient next-nearest-neighbor hopping t' < . 
The double-layer structure of YBCO or BSCCO is not important for the resonance to form. The 
resonance energy u res and spectral weight at optimal doping come out comparable to experiment. 
The observed qualitative behavior of uj res with hole filling is reproduced in the underdoped as well 
as overdoped regime. A second, much broader peak becomes visible in the magnetic excitation 
spectrum if the 2D wave- vector is integrated over. It is caused by excitations across the maximum 
gap, and in contrast to the resonance its energy is almost independent of doping. At energies above 
or below tj res the commensurate resonance splits into incommensurate peaks, located off (tt, tt) . 
Below ujres the intensity pattern is of "parallel" type and the dispersion relation of incommensurate 
peaks has a negative curvature. This is in accordance with recent INS experiments on YBCO . 



PACS numbers: 71.10.Fd, 74.25.Ha, 74.72.Bk, 75.20.Hr 



I. INTRODUCTION 

The study of magnetic excitations plays an important 
role in the ongoing attempt to understand the physics of 
High- Temperature Superconductors. A key observation 
is the so-called "41 meV resonance" from inelastic neu- 
tron scattering (INS) experiments phm ■ In supercon- 
ducting optimally doped YBa2Cu306+j, (YBCOq+ v ) a 
sharp peak occurs in the magnetic structure factor at the 
antifcrromagnetic (AF) wave vector q = (tt, tt) and en- 
ergy 41 meV . It appears resolution limited in energy and 
therefore is described as an (undamped) <5-peak. This is 
not expected for a d-wave superconductor, since the den- 
sity of states is finite and the resonance energy u> res rs 
40 meV is not small compared to 2A° , with the maxi- 
mum gap A ~ 30 ... 40 meV J5| . When temperature is 
raised through T c m 93 K into the normal state, the reso- 
nance vanishes. The main effect of underdoping J6|-|9| on 
the resonance in the superconducting state is a continu- 
ous reduction of its energy, as far as uj res ss 24 meV for 
the most underdoped samples T c rj 50 K . The resonance 
also gains spectral weight with underdoping. In contrast 
to the optimally doped case it persists into the normal 
state above T c , where the pseudo-gap regime is found. 
Recently a resonance has also been observed pJ-O] in 
Bi 2 Sr 2 CaCu20 8+( 5 (BSCC0 8+(5 ) . Its energy 43meV in 
the optimally doped sample is comparable to the case of 
YBCO . If experiments on YBCO and BSCCO are put 
together, cu res seems to follow T c , i.e., it is maximal for 
optimal doping and is reduced in under- as well as over- 
doped compounds llpT) • 

In this paper we report theoretical calculations of the 



magnetic excitation spectrum in YBCO and BSCCO . 
Our starting point is the doped Mott insulator, described 
by the t- J-model. We follow the standard procedure 
of introducing auxiliary "slave" particles and treating 
the problem in mean-field theory. In the resulting ef- 
fective theory the d-wave superconducting phase and 
the pseudo-gap regime of underdoped systems are rep- 
resented as spin-charge separated states. The dynamical 
magnetic susceptibility is obtained from a Fermi liquid 
of strongly renormalized quasi particles that carry the 
spin. The "41 meV resonance" is interpreted as a collec- 
tive spin-1 excitation, it arises from a particle-hole (ph) 
bound state of these quasi particles. We find results in 
good qualitative agreement with the neutron-scattering 
experiments. In particular, the behavior of the resonance 
energy w res with hole filling is reproduced for underdoped 
as well as overdoped systems, and we obtain reasonable 
absolute values for u) res and spectral weight of the reso- 
nance at optimal doping. Our findings are discussed in 
detail from the doping dependent band structure of the 
quasi particles. 

The concept of the resonance coming from a ph- 
bound state has been put forward in several studies us- 
ing "slave" -particle schemes for t-J and Hubbard models 
]13| fl7| , a Hubbard-operator technique |L8| , approaches 
based on BCS theory [fpl-E3| , and self-consistent treat- 
ments of spin fluctuations in the Hubbard model (FLEX) 
|4j-§§ or spin-fermion model [|7|J28) . In the SO (5) ap- 
proach py,p(| , on the other hand, the resonance is a re- 
sult of a bound state in the spin-triplet particle-particle 
(pp) channel, which couples to the magnetic ph-channel 
in the superconducting state. In Ref. |H| we studied 



the contribution from the pp-channel within the present 
slave-particle scheme and concluded that it cannot give 
rise to a resonance below 2A° unless unreasonable pa- 
rameters are used. A similar conclusion has been given 
in Ref. || . 

The resonance is connected to incommensurate struc- 
tures in wave-vector space. Above and below the res- 
onance energy uj res a splitting of the single peak at 
q = (71", 7r) into four peaks slightly displaced from (tt, tt) is 
observed 33-381 . When energy is raised from u) res these 
follow a dispersion w(q) similar to AF spin waves; how- 
ever, the peaks are very broad. Below ui res four well sepa- 
rated peaks are visible, which move away from (tt, tt) with 
decreasing energy and hence are described by an 'upside- 
down' dispersion. Interestingly these peaks are displaced 
from the AF wave vector in direction of the (tt, 0) or 
(0, 7r) points, rotated from the nodal directions by 45° . 
This is the same "parallel" type of incommensurability 
as is known from the L^-zSr^CuC^+y (LSCO) family of 
compounds |3^|I(J , where it has been brought into con- 
nection to the so-called "stripes" . In the present work we 
do not consider the possibility of a combined ordering of 
spin and charge into quasi one-dimensional (stripe-like) 
structures. Nevertheless, below the resonance energy we 
obtain an incommensurate pattern of parallel type. This 
is due to a "dynamic nesting" mechanism in the super- 
conducting state that enhances the intensity at these par- 
ticular points in wave-vector space. The dispersion rela- 
tions of incommensurate peaks are traced back to two 
particle-hole excitation thresholds that vary differently 
with wave-vector. 

Recently the magnetic response has also been stud- 
ied by averaging the neutron-scattering data over the 
2D Brillouin zone [p p^j4l| ] . Besides the resonance, the 
resulting local magnetic excitation spectrum Imx2£> (w) 
shows a second, broad feature at an energy above the res- 
onance, which depends only weakly on the doping level. 
Within our calculation this feature is naturally explained 
from particle-hole excitations across the maximum d- 
wave gap A . Their energy QJhump ^ 2A° comes out 
almost independent of doping. 

The paper is organized as follows: In Sections [F| to 
IV the mean-field theory for the t— t'-J-model is derived 
and some basic implications are reviewed. The magnetic 
resonance at the AF wave vector (71", tt) is considered in 
Section [V] for a single CuC>2 plane. Section [Vj presents 
results for the magnetic response in wave-vector space. 
We consider the crossover from commensurate to incom- 
mensurate response and the dispersion of incommensu- 
rate neutron peaks. In Section VII wc take into account 
that YBCO and BSCCO are actually bi-layer materials 
with two coupled CuC>2 planes per unit cell. The split- 
ting of the susceptibility into two modes is calculated. 
The single-layer model considered in the previous sec- 
tions serves as an effective model for the odd ( "acoustic" ) 
mode, where the resonance is observed. In this section 



we also discuss the above-mentioned local susceptibility 
ImX2Z)(w) . A summary is given in Section VIII . 



Some of the results have been presented briefly in Refs. 
[ pl| , p2[ . Work of other authors is further referenced in the 
respective sections. 



II. MODEL AND MEAN-FIELD THEORY 

We study the t- J-model on a simple square lattice of 
Cu-3d orbitals for each of the two CuC>2 layers in YBCO 
or BSCCO : 
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(1) 



In the subspace with no doubly occupied orbitals, the 
electron operator on a Cu-lattice site v is denoted cw 
with spin index a = ±1 ; S v is the spin-density oper- 
ator. A Cu-site is specified through v = [i,l], where 
i = l... Nl indicates the Cu-position within one Cu02- 
plane and I — 1,2 selects the layer in the bi-layer sand- 
wich. t vv i denotes the effective intra- and inter-layer Cu- 
Cu-hopping matrix elements, and J vv i the antiferromag- 
netic super exchange . 

To deal with the constraint of no double occupancy, 
the standard auxiliary-particle formulation is used, 



btU 



(2) 



The fermion /J CT creates a singly occupied site (with spin 
a), the "slave" boson 6J, an empty one out of the (un- 
physical) vacuum b v \0) = fv<r\0) = 0. The constraint 
now takes the form 



h=b%+y,ft<rU=i 



(3) 



Using (0), local operators can be written in fermions, in 
particular the particle and spin density read 



iptipu , S v = -$ItiI) v 



(4) 



Here spinors 



^=te)^=(j; 



have been intro- 



duced, with Pauli matrices f and h = 1 . 

In order to derive a mean-field theory the constraint 
(||) is relaxed to its thermal average (Q v ) = 1 . Together 
with the number x of doped holes per Cu-site, it fixes the 
fermion and boson densities to 



(1 - x) = (1>M , x = (blK) 



(5) 



These are adjusted by chemical potentials fi b ,^ ■ Using 
a coherent-state path integral the partition function is 
now represented by the action 



S = S° + S t + S J + s h 



(6) 



with 



S°= / dTJ2{b»(8 T -H b )bv + Mb- -H f )*!>»} 
ft 3 

Jo , 
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S J = / (ir-Vj w -5A- , S h = - drViS, 
Jo l „ „, Jo 



A magnetic source- field ft = /i„(t) has been added here, 

0=l/k B T. 

A mean-field decomposition of the interaction terms 
S f , S J is achieved via Feynman's variational principle |43J ] 
for the free energy _F , 



f3F<<£>[S] , <£>[S] = (S-S)- In Z 



(7) 



The effective action 

S determines Z = fD[ip, tp,b, b] exp(— S) and thermal 
averages 

(d) = i [r>ty,ip,b,b]exp(-s)d 



For S we make the quadratic ansatz 
S = S°+ I 'dld2J^ 1 T{ 2 ^ 2 

+ hTi 2 b 2 + (^iAi 2 ^2 + h.c. 
dlrhiSi 



(8) 



with a shorthand notation 1 = {v\,Ti) , Jdl = 

Si/i /o^ ri collecting site v and time r indices. T{ 2 and 
A12 are matrices in spin space, e.g., A i2 = A^ 2 = 
A1\Zl ( T ii r 2) ■ ^ consists of quadratic terms for fermions 
and bosons, which represent all possible mean-field de- 
couplings of the interactions 5^ and S J . 

The expectation value (5 — S) in Eq. (0) is calculated 
using Wick's theorem, and from the vanishing variation 
S^f[S] — we obtain equations for the self-consistent 
parameters, 



3* = -tia^^i) 

T^ " = - ti 2 (&2&l)<5a 



A^ = - Jxi\Y l {{r^Y {r^Y) 



mi = hi - I d2 J^Sz) 



(9a) 
(9b) 



(9c) 
(9d) 



r M denotes a Pauli matrix, rj its transpose. The effec- 
tive hopping T b of bosons as well as the first contribution 
to the hopping T* of fermions stem from the decoupling 
of (£*) in Eq. @j . The Heisenberg term (S J ) is factor- 
ized through Wick's theorem into contributions to the 
local magnetic field m , the fermion's hopping T* and 
pairing amplitude A . These correspond to analyzing S J 
in the direct particle-hole (ph) channel of fermions, the 
exchange ph, and the particle-particle channel, respec- 
tively |44J . Further below, when the source field ft is set 
to zero, these will be restricted to a Resonating Valence- 
Bond (RVB) amplitude Tl v , ~ (Jv\$v'\) and spin-singlet 
pairing A vv > ~ (f^fu'i) ■ 

The approximate free energy is the functional $ at 
its stationary point, /3F appr = \&[S stat ] , with the action 
gstat gj ven j-jy gq (|h an( j gq S> (|g|) _ phe dynamical mag- 
netic susceptibility then follows with ^[S' 5 *' 2 *] /<5fti = 
-(&) as 
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Using Eq. (9d) we obtain the usual 'RPA-like' expression; 

(10) 



in a matrix notation it reads 
X = [1 + Jx"' 



X 



The irreducible part is identified as 



A12 



S(S 2 ) 
Srrii 



= (S1S2 




(11) 
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The 1st term on the r.h.s. comes from the rn which ap- 
pears explicitly in (§2) through S , Eq. (N) ■ For a vanish- 
ing source field ft, = it is given by the unrenormalized 
fermion bubble contained in Fig. m. Since the operator 
S involves only fermions, no boson excitation occur in 
X lvr ■ The 2nd term in Eq. (fyl) stands for all contribu- 
tions from the implicit m-dependence of (S?) through the 
other mean-field parameters T b ,T^,A. Using Eqs. (J9a|- 
|9q ) , it gives rise to vertex corrections. These are shown 
in Fig. [l] (bottom). 

In the following we set ft = and consider paramag- 
netic phases m — , which are symmetric with respect to 
lattice translations within a CuC>2-layer and exchange of 
the layers. In going to wave-vector space, the site index 
v = [i,l], with in-plane site i = 1 . . . Nl and layer index 
I = 1, 2 is replaced by the wave- vector p = (k,p z ) . That 



is, 



fila 



y/2Ni 



■£■ 



1 * ^ ™ x * a; "T" ™ J/ *^y / 



E 



./per 



and similar for boson operators. Here k runs over the 
usual 2D Brillouin zone, and p z = 0, 7r corresponds to 
even, odd linear combination of layer orbitals. 





with 



FIG. 1. Top: Irreducible part x l " f° r vanishing source 
field h — . Full lines are fermion Green's functions, bo- 
son excitations do not enter x' rr a t mean-field level. Bot- 
tom: Bethe-Salpeter equation for the vertex function in \ %rr , 
dashed lines stand for the Heisenberg interaction J . The bare 
vertex (i.e., the unrenormalized bubble) represents {SiS2) conn 
in Eq. fill) . 



(12) 



The exchange interaction is decomposed as 



1 , = I 



w 



SwJ5<i,j> + (1 - 5w)J 5i 



It consists of an intra-layer component J for nearest 
neighbors < i,j > and an inter-layer coupling J for 
i = j . Thus the pairing A of fermions, Eq. (pa), involves 
an intra-layer part, which we restrict to singlet pairing 
with d-wave symmetry and equal amplitude and phase 
in both layers: (futfu+xl) = -(flitfu+yl) > l = l > 2 - 
A contains as well an inter-layer amplitude (/li-f/btj.) • 
In p-space Eq. (Bq) then becomes A\% — > A aa (p) = 
crS-a' o-Ap , with the gap function 



A, 



A c 



cos 



\{k x ) - cos{k y )\ + A-^e*" 



and the maximum in-plane and inter-plane gap 
A° = |JA,A X0 = |J X A X , 

A = (/iif/ii+44.) — (fnihi+xy 

A = (flilf2n) ~ (flilf2i-[) 



(13) 



(14a) 
(14b) 



As will be explained at the end of this section, in the 
interesting range of temperature and hole filling x the 
bosons may be treated as almost condensed. That is, 
the hopping rate {b v b v i) rs {b v b v ) — x is independent 
of v, v' and given by the hole density x via Eq. (pf) . 
The first term in the fermion hopping Eq. (BH) becomes 
— t vv '(b v b U ') — ► —xt vv i . It describes the propagation of 
fermions with the small probability x of finding an empty 



site. The second term in Eq. (9b) involves induced hop- 
ping amplitudes on nearest-neighbor bonds, which we as- 
sume equal in amplitude and phase on each bond within 
a layer (uniform RVB), {fu-\.fii +x -\) = (fmfu+yi) , 
I = 1,2. The fermion hopping Eq. (9b) now turns into 

Tt°° \v) = 5 aa ,Ti , 



rpfcra' 

± \1 



T/ = - 2i[cos(fc x ) + cos(fc y )] 

— At cos(k x ) cos(ky) — t (k)e 



(15) 



and 



t = x t A — Jy , t = xt 

8 

^(k)=.T^(k) + |j^x ± 



X — (/lJl/lj+xl) + (flilfli+xl) 

X =(/l!t/2'»t) + (/l'»l/2il) 



(16) 



(17a) 
(17b) 



For the bare hopping elements we assumed a nearest and 
next-nearest neighbor overlap t and t' within a layer, 
and an inter-layer S,p6[ hopping i x (k) = 2t ± [cos(k x ) — 
cos(fc y )] 2 + ig" . For bosons the effective hopping is de- 
rived similarly, T\ 2 — » E p . The result is given at the end 
of this section. 

The mean-field Hamiltonian Eq. (ph now reads 



S = f dr J2 \b p (d r + n p )b p + J2 



fpa(d T +£ p )f pa (18) 



+ [ A p(/pt/-p1 - Ipll-p]) + h-c. 



It consists of free bosons with dispersion £l p = J*—// 6 and 
BCS-fermions with e p — T* — \J and gap function Eq. 
( |l3| ) . After Bogoliubov transformation we obtain 'quasi- 

fermion' energies E p = Je^ + A 2 , and the mean-field 

equations (fl4|), ( |17| ) become 

Nt 2—< 
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E, 



f tanh(/3S p /2) (19a) 
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jL ^^tanh(^ p /2) 



2A^ L ^ B. 



(19c) 



with phase factors 7(k) = cos(fc x ) + cos(k y ) , y(k) = 
cos(fca;) — cos(fcjj) . The last equation is the particle num- 
ber constraint (|5|) . 

The magnetic susceptibility Eq. (|l^) is isotropic in 
spin-space for h = fn = , and takes the usual form 



XpH 



4 rr M 



l + jpxjrM 

where J p is obtained from Eq. (|l2|) as 



J p = 2 J [cos^) + cos(g y )] + e lp * J J 



(20) 



(21) 



In experiment the magnetic response ^ meas is measured 
as a function of the wave vector (q, q z ) , which spans the 
3D Brillouin zone of the bi-layer material. It is given by 



X meas (q,q Zl Lu) 



(22) 



0.5^ 



(ffMs) 



2fd \ / m • 2/d 



Xp( w )L=o cos (o?*)+ X P (w)| w sin (-ft.) 



with p = (q,p z ). d denotes the spacing of layers in 
the double-layer. The even (p z = 0) and odd (p z = it) 
mode susceptibilities correspond to the in-phase and anti- 
phase combination of spin fluctuations in the planes. For 
the irreducible part Eq. ( |ll[) we take the bare bubble 
x£J/(t,t') = (S*(T)S*,(T')} conn . The vertex corrections 
depicted in the bottom of Fi g, [l] can be savely ignored. 
As we have discussed in Ref. |3l| they have no significant 
effect in the interesting energy range < u> < 2A° . With 
the effective Hamiltonian Eq. ( |18| ) we get the expression 
known from BCS theory 
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(23) 
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f(s'E') - fjsE) 
+ sE - s'E 1 + ?:o H 



Here e = Ep , e' = £p+ P , and similar for A , E . f denotes 
the Fermi function. 

We close this section with a remark on Bose conden- 
sation. From Eq. ( pa]) the boson dispersion is fl p = 
— 2tx[cos(k x ) + cos(ky)] — /i b , where for simplicity t' — 
t 1 - = . Near the band minimum k — this becomes 
Cl p w fl + k 2 J 2m b , with the mass 1 / m b — 2t\ ■ From 
the solution of Eqs. (19) we get values around x ~ 0.4 , 
i.e., m w 1/t. In two dimensions free bosons do not 
condense at finite temperature T > , however, the cor- 
relation length of the propagator (b v b v >) grows exponen- 
tially for T below T% E = 2irx / m b w 2irxt . In the x 
and T range we are interested in, T <C Tg E , and the 



bosons can be considered almost condensed, i.e., ft 
and (b v b v i) ss (b^b^) = x for any z/, v' . 







III. PHASE DIAGRAM 

The slave-boson mean-field theory has been put for- 
ward in numerous papers [}48|-p2| , originating in the Res- 
onating Valence-Bond (RVB) idea |p3[ . In this section 
we review the phase diagram and briefly discuss some ex- 
perimental implications in the superconducting phase at 
T — ► . For simplicity a single CuC>2-layer is considered, 
with t = 2 J . Fig. shows the phase diagram, derived 
from the numerical solution of Eqs. (n9|) . It resembles 
those given in the literature plpq] , except at very small 
doping x — ► , where our assumption of almost con- 
densed bosons becomes incorrect. We also ignored the 
staggered- flux phase reported J56fl for small x . These 
simplifications do not affect our results for the magnetic 
excitations. 

Ignoring the lines labeled 'AF' for the moment, the 
phase diagram shows two transition temperatures Td 




0.2 0.25 

Hole filling x 
FIG. 2. Main figure: Mean-field phase diagram for a 
single Cu02-layer with t = 2,1 . Td and Tbe denote transition 
temperatures to d-wave pairing of fermions and condensation 
of bosons, respectively. T c = min{Td,TB_B} is the physical 
T c for bulk d-wave superconductivity. The lines labeled AF 
indicate the magnetic instability of the bare (dotted line) and 
renormalized theory (continuous), see text, Sec. nV\. Inset: 



Maximum gap A 
T = 0. 



of fermions as function of hole filling at 



and Tbe ■ At temperatures T > T^ the fermions move 
in the CuGvplane with an effective dispersion e(k) = 
— 2t[cos(k x ) + cos(ky)] — fif , where \ m t — %t + |Jx 
is finite (uniform RVB phase). At T = Td they undergo 
a transition to the d-wave paired state with order pa- 
rameter A(k) = ^-[cos(k x ) — cos{ky)\ . This spin-gap 
phase is characterized by a gap A(k) for spin excita- 
tions, and is associated with the spin-gap (or pseudo- 
gap) regime observed in the normal state of underdoped 
cuprates. The bosons, on the other hand, show Bose 
condensation (b u ) ^ at T = Tbe ■ Bulk supercon- 
ductivity is present for T < T c = min-fT^, Tbe} , where 
(cv\c v >i) = (&£)(&£') (futfu'i) + 0. In two dimensions 
Tbe is identified with Tg B = 2-ixxt , which yields a very 
large slope of the Tbe line in the phase diagram. Fluc- 
tuations of gauge fields around the mean-field solution 
are expected to reduce T BE to reasonable values |5J| . 
The Tbe line sketched in Fig. g corresponds to that sit- 
uation, with a maximum T c at an optimal doping value 
x op t ~ 0.15. The similarity of spin-gap and supercon- 
ducting (SC) phase in mean-field theory naturally ex- 
plains the persistence of the magnetic resonance into the 
spin-gap regime, although line-shape (damping) effects 
are missing. 

In the following we focus on the SC state at T — ► , 
which is reasonably reproduced by mean-field theory: In 
cuprate superconductors the underdoped region x < x opt 
shows unusual behavior of T c , the superfluid density, and 
the maximum gap as function of hole filling. The super- 
conducting T c increases with doping, T c — Tbe ~ x . 



The superfluid density p s = n s /m is given by the con- 
densate density p b = n b /m b ~ x of bosons for small 
x , according to the Ioffe-Larkin formula (57| . Thus 
p s ~ x , and the well known experimental observation 
p8fl p s ~ T c follows naturally. Tunneling and photoe- 
mission spectra are described by the Green's function 
Gm'(t,t') — (T t c u ]-(t)cI,Jt')} , where c is expressed by 
Eq. (0) . In the superconducting state, where bosons are 
condensed in k = , G splits into a coherent and inco- 
herent part, G = xG ferm + G mcoh , where G ferm is the 
propagator of fermions. Thus the superconducting gap is 
given by the d-wave pairing gap of the fermions. Its dop- 
ing dependence at T = is shown in the inset of Fig. g. 
When x is reduced from x op t the maximum gap actually 
increases (whereas T c decreases) , as is seen in experi- 
ment |59|-|62[| . At optimal doping mean-field theory gives 
A « 0.3 J rs 40meV, which compares reasonably with 
experimental values || . On the overdoped side x > x opt , 
T c = T d , and we get the BCS-like result A - T c . 

Recently an alternative slave-boson formulation has 
been proposed (6^Q] , which extends the SU(2) symme- 
try in particle-hole space of the 1/2-filled model |35| to 
the hole-doped case. Within mean-field theory the super- 
conducting state at T — ► appears to be similar to the 
more conventional U(l) formulation we are using here, in 
particular is the magnetic spin-response the same. 



IV. EFFECTIVE INTERACTION 



Jxaf ■ Coming from high temperature or doping, an AF 
instability is indicated by (xaf) -1 —> 0. xaf has been 
calculated numerically, the resulting phase boundary is 
shown in Fig. g as a dotted line labeled 'AF'. Our result 
is similar to the phase diagram obtained in Ref. ]69| . 

Apparently, at zero temperature, AF order occurs at a 
quite high hole concentration x° c ss 0.22 , which is totally 
inconsistent with experiment (x c sw 0.02). Furthermore, 
the study of magnetic properties in the paramagnetic 
phase is bound to x > x° c , i.e., the overdoped region. 
The high x° c is an artifact of the mean-field approxima- 
tion. Within the gauge-field approach it has been shown 
Q that the AF ordered state at 1/2-filling x = is 
quickly removed for x > . Furthermore it is known 
p3|,[71 72 1 that the interaction of spin- waves with doped 
holes destroys AF order at a small x c . In order to treat 
underdoped systems we include these physics in a phe- 
nomenological fashion. We assume a renormalization of 
the magnetic interaction J p in Eq. fl20|) , such that x° c is 
reduced to some x c ~ 0.03 . The model is to replace 



J — > aJ , a = 0.35 <-> x r 



0.03 



(26) 



in Eq. (^TJ) , J 1 - stays unchanged. The actual value of 
a is equivalent to choosing a specific critical doping x c . 
As long as x c is physically reasonable (< 0.05), results 
do not depend significantly on x c . Using J — v a J in Eq. 
(E4T) we get a new AF phase boundary, which is shown 
in Fig. as the continuous line 'AF'. 



So far we have not considered the possibility of anti- 
ferromagnetic (AF) order. It is known that wide areas of 
the mean-field phase diagram are unstable to AF order 
Ipq p9| . We determine the AF phase boundary from the 
diverging correlation length £af , which is extracted from 
the static (u> — 0) susceptibility. For a single Cu02-layer 
Eqs. ©-© reduce to 

X[q ' UJ ' l + 2J[cos(q x )+cos(q y )}x irr (ci^) [ ' 



with x m units 47 of (ff/is) 2 ■ The irreducible part 
X 4rr (q, w) = Xp r {u) is given by Eq. (||) withp = (q, p z ) , 
p = (k,p z ) and arbitrary p z . In Eq. (p3]), the in- 
ternal summation is now over the 2D Brillouin zone, 
2l^Ep -> l^Ek = X! 7r d 2 fe/(27r) 2 , and dispersion 
and gap function become 



e = e(k) = 
A = A(k) 
E = E{k) 



Ti 



J 



(i- with t = ; e = e(k + q) 
Ap with A ±0 = ; A' = A(k + q) (25) 
^e(k) 2 + A(k) 2 ; E' 



E(k- 



A p , T/ have been given in Eqs. (|Tjjrj,|TJ) . At the Neel 
wave- vector Q = (n, n) the static susceptibility takes the 
value xaf — x(Q?0), and for wave vectors q close to 
Q we get J X (q, 0) = 1 / [(q - Q) 2 + ££] with & F = 
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FIG. 3. Main figure: Antiferromagnetic (AF) corre- 
lation length (af at T = in units of the lattice spacing. 
Continuous line: £,af{x) for an effective interaction a — 0.35 
(see text). The vertical dotted line indicates the AF insta- 
bility at x c — 0.032 . Dashed line: Function 0.2/^/a; — x c , 
fitted to £af(x) for 0.033 < x < 0.1 . Inset: Static suscepti- 
bility X^iO) (i n arbitrary units) as a wave- vector scan over 
(7r,7r) in the 2D Brillouin zone, for x — 0.033 slightly above 
x c . Shown is a scan in (-7r,0)-direction, scans in other direc- 
tions look similar, indicating a commensurate AF transition. 
a — 0.35 and T = as above. 



To give another argument in favor of the simple inter- 
action model, especially the a being independent of dop- 
ing, we consider the correlation length £af (x) as function 
of hole filling. Fig. || shows £af > calculated with J — v 
a J . It diverges at a x c — 0.032 and decreases rapidly 
with additional doping, following £,af{x) ~ 0.2/^/a; — x c . 
This behavior is consistent with neutron-scattering mea- 
surements p3[ on LSCO and results from high tempera- 
ture series for the t- J- model Q . The function ~ 1/v 7 ^ 
represents the average distance of doped holes and has 
been used in [J73| to interpret the data. Finally it is noted 
that the AF transition at x c occurs at the Neel wave vec- 
tor Q = (71", 7r) , i.e., is commensurate. This is shown in 
the inset of Fig. ||. 



V. THE MAGNETIC RESONANCE 

This section presents results for the magnetic response 
at the antiferromagnetic (AF) wave vector Q = (w, ir) . 
We consider a single Cu02-layer with a nearest and next- 
nearest neighbor hopping t = 2J and t' — — 0.45£ , appro- 
priate ||75| , [46|l for YBCO . Effects specific to the bi-layer 
structure of YBCO and BSSCO, namely the splitting of 
the magnetic response into acoustic a nd o ptical modes 
will be discussed separately in Section VII . 
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FIG. 4. Magnetic response at the AF wave vector (71", n) 
for a single Cu02-layer near optimal doping, x = 0.12. Pa- 
rameters are t — 2J , t' = — 0.45£ . Main figure: Imag- 
inary part of the bubble Y^ rr m Eq. (J2J). Shown is the 
superconductor at T = (cont. line) and the normal state 
T = 0.2J > T c (dashed line). The dashed-dotted line is the 
inverse Stoner-factor K (see text) for T — , scaled x (—5) . 
Inset: Imaginary part of the resulting susceptibility x, Eq. 
(||), for T = (cont. line) and T > T c (dashed line). The 
sharp peak visible for T = is actually a <5-function, broad- 
ened by a small damping used in the numerical calculation. 



A. Results for the AF wave vector (n,n) 



The dynamical susceptibility is obtained from Eq. (|24| ) 
with the effective interaction J — > a J , a — 0.35 intro- 
duced in the preceding section. The integration in Eq. 
( P3| ) is performed numerically on a 5000 x 5000 lattice in 
k-space, with the infinitesimal i0 + — > i2T replaced by a 
small finite damping 2r = 0.001 J. 

Calculations for a fixed hole filling x = 0.12 near 
optimal doping, based on mean-field parameters from 
the self-consistent solution of Eqs. (19) are shown in 
Fig. [|. In the superconducting (SC) state at T = 
the imaginary part Imx irr (Q, w) of the irreducible bub- 
ble is characterized by a gap up to a threshold en- 
ergy J7o rs 0.54 J , with a step-like van Hove singularity 
(v.H.s.) at the onset of spectral weight at w = Slo ■ A 
peak at to s=s 2A° = 0.72J is remnant of the density of 
states of the d-wave superconductor. The correspond- 
ing real part is shown in Fig. as the inverse Stoner- 
cnhancement factor K(Q,w) = [1 — a4JRex irr (Q,^)] ■ 
By virtue of the Kramers-Kroenig transformation, the 
step at the threshold f2rj in Imx irr turns into a log- 
singularity in Kex trr , and K(Q,uj) crosses zero at an 
energy uj res = 0.5 J < ^o within the gap. This leads to 
an undamped (5-like resonance at uj res in the magnetic 
response Imx(Q,w) , as is shown in the inset of Fig. 0. 
The position ui res — 0.5 J w 60meV is not too far off 
the « 40meV observed in optimally doped YBCO and 
BSCCO . The situation changes drastically in going to 
the normal state T > T c w Td . As is seen from the 
dashed line in Fig. |j (main figure), Imx irr (Q,w) looses 
its structure, in particular the gap vanishes. The cor- 
responding K (not shown in the figure) becomes equally 
structure-less without any zero crossing, leading to a van- 
ishing of the resonance in Im^ in the normal state (inset 
of Fig. |). 

The effect of undcrdoping is demonstrated in Fig. |5|. 
For comparison with experiment Imx in the inset has 
been broadened to an experimental resolution (FWHM) 
of 4r ss 5 meV through a damping i0+ — > 2r = 0.02 J in 
Eq. (|23| ) ■ When x is reduced from x = 0.12 the gap Q.q 
in lva.x irr decreases monotonously, and with it the reso- 



nance at tores ^S ^o nioves to lower energies. u> res reaches 
zero at the AF transition, which occurs at x c = 0.023 
for the parameters used here (x c depends slightly on t' , 
since a = 0.35 is held fixed). Note that uj res moves oppo- 
site to the maximum gap A . The latter increases when 
x decreases (see Fig. S) and is reflected in the peak in 
Imx irr at higher energies w = 2A° = 0.72 . . . 0.80 in Fig. 
H. The spectral weight W — J dwlmx increases |7q] 
when x is reduced, since the system is shifted closer to 
the magnetic instability. For a quantitative comparison 
of W we follow a procedure applied to experimental INS 
data in Ref. p|M| . For x = 0.12 the flat normal-state 
spectrum is subtracted as a background from the T = 
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FIG. 5. Magnetic response at wave vector (ir, 7r) as in 
Fig. ||, for optimal to underdoped hole filling in the super- 
conducting state T = . Main figure: Imaginary part of 
X lrr ■ Curves are identified by the respective onset of spectral 
weight (threshold), which is 0.54 J for x — 0.12 and decreases 
for x = 0.08, 0.06, 0.04, 0.023 down to 0.09J . Inset: Imag. 
part of the susceptibility x ■ Here the peaks are broadened to 
an experimental resolution (FWHM) of 4F ~ 5 meV . 



curve shown in Fig. 0, inset. Integrating only the pos- 
itive part of the resulting difference spectrum gives the 
weight AW of the resonance compared to the normal 
state. We find AW = 1.55 /x^, which agrees well with 
optimally doped YBCO . With reducing x the resonance 
also develops some intrinsic damping. The step-height 
at the threshold f2o m Imx 4rr , which is responsible for 
the (5-like resonance, decreases and eventually vanishes 
around x = 0.09 . 
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FIG. 6. Magnetic response at q = (n, n) for optimal to 
overdoping at T = . Parameters as in Fig. H , except hole 
filling. Main figure: Bubble spectrum Imx lrr for x — 0.12 
(curve a), 0.14 (b), 0.18 (c), 0.24 (d), 0.30 (e) . Inset: Cor- 
responding Imx , broadened to 4F w 5 meV . 



The effect of overdoping is presented in Fig. M. When 
x is increased from 0.12 up to 0.3, the ph-threshold fi 
in Imx irr and with it the position of the resonance first 
grows, but ui res starts to decrease around x = 0.18 . This 
trend is consistent with recent INS experiments [jTlJ on 
overdoped BSCCO, where the resonance appeared at an 
energy uj res reduced from the optimally doped case JlfJ . 
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FIG. 7. Left panel: Particle-hole excitation threshold 
fio and resonance position uj res as function of hole filling in 
the SC state, ilo is given by Eq. (127]) , and uj res is the energy 
of the peak maximum in Imx , extracted from plots like Figs. 
pL H (insets). Right panel: Comparison of Qq , the chemical 
potential jj, — —\jjr\ of fermions , and the maximum gap A . 



Fig. M (left) summarizes the doping dependence of the 
resonance. u> res is always located slightly below the 
threshold f^o to the damping particle hole continuum. 
Near the magnetic instability |77|] at x = x c = 0.023 and 
around optimal doping x ~ 0.1 . . . 0.2 , the resonance is 
a ^-function, well separated from the continuum. The 
x-dependent threshold is given by 



2|^ / | v /cr(2-cr) for a < 1 
fio = < (27) 



!/-<■ 



/i 



for a > 1 



with a = (A ) 2 / 8xt'n f . For x < x w 0.09 it is a > 1 
and therefore fin = 2|/i' | . That is, in the underdoped 
regime the resonance energy follows the chemical poten- 
tial I/** | = —\v of the fermions, and thus increases with 
hole filling. The gap A , on the contrary, decreases. 
This is illustrated in the right panel of Fig. [7|. Around 
x = 0.13, where \\J\ = A , a crossover into the over- 
doped regime occurs, where $7o « 2A° for large x . The 
increase of tv res turns into a decrease. Compared to ex- 
periment |ll| the latter is too weak; this is due to A 
(and T c = Td) decreasing too slowly with x in the self- 
consistent mean-field calculation. 

It has been pointed out above that the resonance de- 
pends on the d-wave gap to be present, i.e., the supercon- 
ducting or spin-gap state, and a sufficient next-nearest 
neighbor hopping t' < . The influence of t' becomes 
apparent, if the bubble x* rr and the resulting magnetic 
spectrum Imx at T = are re-calculated for t' = (not 
shown in the figures). The response Imx no longer con- 
tains a resonance-like peak, although the spectral weight 
is conserved. Only for small x — > x c at the magnetic 



instability Imx develops a Goldstone (Bragg) peak. The 
bubble spectrum at the threshold Slo = 2|/i^| now follows 
Imx* rr (Q,a;) ~ y/oj — Q Q(ui — fi ) for the whole range 
of x , and a step-like v.H.s. never appears. 



B. Discussion 



The results presented above compare well to neutron- 
scattering experiments ||pf : The variation of the res- 
onance energy oj res with hole filling x is reproduced in 
the underdoped [||j7|,|l| and the overdoped @,0,0 
regime. At optimal doping the resonance appears reso- 
lution limited (as a delta function) in the SC state only 
[g,[| , and its energy and spectral weight are compara- 
ble to the values measured in experiment |8j . In the 
underdoped regime it is obtained also in the spin-gap 
phase [pp,|l2]| above T c ; the observed line-shape (damp- 
ing) is not reproduced in mean-field theory. YBCO and 
BSCCO are bi-layer materials, i.e., consist of two coupled 
CuC>2 planes per unit cell. However, for the resonance to 
emerge the bi-layer structure is not important. Rather, 
in optimally and slightly underdoped systems it depends 
on a hole-like Fermi surface (i.e., a sufficient t' < 0) and 
a finite d-wave (spin-) gap. As will become clear in Sec- 
tion 



VII this conclusion is not altered if the double-layer 



structure is taken into account. 

The slave-particle approach reproduces quite satis- 
factorily the resonance energy u> res in the underdoped 
regime, where uj res is not connected to the maximum 
gap A , see Fig. [7|. u res is given by the pole of Eq. 
(pi|), i.e., the energy of the bound state in the particle- 
hole channel of the fermions. In the underdoped regime 
^res ^ 2|/x^| follows the chemical potential \^\ = — yj . 
\i' refers to quasi particles (the fermions) which emerge 
from the mean-field description of the t-J-model and 
arc strongly rcnormalized. They propagate with hopping 
matrix-elements t = xt + | J\ ~ xt + 0.15 J , t 1 — xt 1 , 
and hence Fermi velocity vp ~ xvp , that are reduced 
from the bare parameters by the small Gutzwiller factor 
x < 0.15. The latter mimics the reduced phase space 
due to local correlations in the doped Mott insulator. 
Accordingly |/i^| comes out small enough, such that in 
the underdoped regime \/x^\ < A , and oj res < 2|/u^| 
is determined by \^\, which increases with hole filling 
x . In contrast, if unrenormalized quasi particles (QP) 
are assumed with bare t,t',vp, the chemical potential 
|/i| ~ -\^\ ^S> A , and uj res < 2 A is connected to the 
gap for almost all x , which decreases with x . To achieve 
a satisfactory result from unrenormalized QP the effec- 
tive (residual) interaction [|78| has to be made doping 
dependent. This is the case in theories based on the 
spin-fermion model [E7U2q| , where the coupling constant 
is controlled by a magnetic correlation length, which can 
be chosen ^-dependent. 
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FIG. 8. Origin of the threshold for particle-hole (ph) 
excitations with wave vector (it, 7r) , corresponding to the 
spectra in Fig. g. Main figure: 2D Brillouin zone (BZ) 
with the underlying normal-state Fermi surface (FS) for 
t — 2 J, t' — — 0.45t as dashed lines. The dashed arrow in- 
dicates a ph-excitation with q = (n,n) at a minimal en- 
ergy fio = . In the superconducting state the FS col- 
lapses to nodes (indicated as dots) , and fio > . The full 
arrow connects two constant-energy lines of quasi-particles 
E(k) — fio/2 . Inset: ph-excitation energies fi(q, k) of 
the superconductor for q = (tv, tt) in the upper right BZ/4 . 
Shown are the minima f2o as dots and the first 5 higher ener- 
gies as lines. Line distance is 0.05 J . 



C. Properties of the ph-threshold fio 

In the remainder of this section we study in some detail 
the origin and qualitative properties of the threshold Qq 
in the bubble spectrum Imx* rr (q, w) at q = (tt, it) ■ The 
presence of this threshold and the step- like onset of spec- 
tral weight at ui = fio in the superconducting state lead 
to the sharp resonance. We consider \mx irr at T = , 
which reads from Eq. (|23|) with (Esl) for uj > 



W rr (q,u;)~^(u;-n(q,k)). 



(28) 



BCS coherence factors have been ignored. It is de- 
termined by the particle-hole (ph) excitation energies 
f2(q, k) = E(k) + E(k + q) of fermions. The pair mo- 
mentum q is set to Q = (tt, it) in the following. Fig. @ 
illustrates the situation for x = 0.12 near optimal dop- 
ing, corresponding to the spectra shown in Fig. H. In the 
normal state T > T c a ph-excitation connects two points 
of the underlying Fermi surface (FS), and the threshold 
fio = minfeil(Q,k) is . In the superconducting state 
a finite f2o is at first sight surprising, since the d-wave 
SC has a finite density of states. However, when the FS 
collapses to nodes, a minimum energy fio > has to be 
paid for ph-excitations with wave- vector Q . It turns out 
that fi(Q,k) is minimal, if E(k) = E(k + Q) = fi /2 , 



with the lines of constant energy E(k) = CIq/2 touch- 
ing the reduced (magnetic) Brillouin zone at 8 points. 
These points are connected by Q as indicated by the full 
arrow in Fig. g. Apparently E(k) — £lo/2 is very flat 
near these points, close to a (dynamic) nesting condition. 
This is due to the bandstructure of the underlying normal 
state, in particular the t' < , and the presence of the 
d-wave gap A(k) . The resulting ph-dispersion Sl(Q,k) 
is shown in the inset of Fig. @. It displays two min- 
ima per 1/4 Brillouin zone with energy Sl(Q,k°) = Hq , 
which fall on the line k® + fc° = tt . The neighbor- 
hood of these minima is quite flat, and Q can be ex- 
panded as Sl(Q,k) « fio + ai(fci) 2 + 02(^2 ) 2 , with fci.2 = 
-j=((k x — fc°) ± {k y — ky)) and relatively small ai,a 2 > . 
Eq. (|28| ) then shows a step-like van Hove singularity 
(v.H.s.) at Q , lm X lrr (Q,uj) ~ @(uj - fl ) / y/aTai ■ The 
value of Qq is given by Eq. ( p7| ) with a < 1 . 

When the hole filling x is reduced, the two minima in 
£7(Q, k) move closer, until they merge at fc° = ky = it/2 
for x = x ss 0.09 , which corresponds to a = 1 in 
Eq. (|27| ) . The step-v.H.s. vanishes on the course (see 
Fig. H)^ since S!(Q,k) becomes increasingly steep in fci 
direction. For x < x we have a > 1 , and Imx irr 
may be approximated by setting t' — , leading to |79[ | 
Imx* rr (Q, u>) ~ v 7 ^ — 51o6(w — Qo) ■ I- e -j the step at the 
threshold Qq has changed into a y/ijj behavior. 



(tt, 0)-direction 



(w, 7r)-direction 



VI. INCOMMENSURATE RESPONSE 



The resonance at uj res = 40 mcV , as well as its rela- 
tive in underdoped samples with reduced uj r es < 40 meV 
is characterized as a single (commensurate) peak at 
q = (tt, tt) in wave-vector space f| f||j • Above the reso- 
nance energy an incommensurate structure has been ob- 
served p3,B6LB7j , with broad maxima following a disper- 
sion similar to spin waves. Recently, inelastic neutron- 
scattering (INS) experiments on underdoped ]9|, p5| , p6| , p0| 
and optimally ||^,^l[] doped YBCO revealed that also 
the magnetic response below the resonance position 
shows incommensuration: Four distinct peaks appear at 
q = (tt ± 6, 7r) and (tt, tt ± 5) , which move away from 
(tt, tt) with decreasing energy, i.e., are described by some 
'upside-down' dispersion J36| , |37[ . The incommensurabil- 
ity is of "parallel" type, peaks are displaced in (tt, 0) 
and (0, 7r)-direction from (tt, tt) , similar to those observed 
||,@ in LSCO . This is not expected in a d-wave BCS 
picture Q,§2| , since at low energies the particle-hole ex- 
citations from node to node should dominate, leading to 
four peaks at q = (tt ± 8', tt ± 5') . In this section we 
demonstrate that a parallel incommensurability actually 
occurs in a range of energies below the resonance. 



30 



=?20 



i 10 

Q. 

E 
< 







1 




• 


1 
1.0 - 

1.1 


" 




___ 


1.2 " 


- 




1 ~~ 


1.3 - 








1.4 


fei£_ 


2*. "-■ 


' l3i 


C&^Z^N 


I 1,1,1, 




-0.2 -0.1 0.1 0.2-0.2 -0.1 0.1 0.2 




-0.1 0.1 
8/271 



0.2-0.2 



-0.1 0.1 

5/2TT 



FIG. 9. Wave-vector q scans of Imx(q, w) at fixed en- 
ergy u, for t — 2J,t' = — 0A5t,x = 0.12 in the su- 
perconducting state at T — . A quasi-particle damping 
r = 0.01J , corresponding to an experimental energy reso- 
lution (FWHM) of 4r « 5 meV has been used, q is mea- 
sured from (tt, tt) 



r): 5 X 



qx, y — tt , 5 — ±y5| + 5% . Left 
column: Scans in (tt, 0)-direction, i.e., S y = . Right col- 
umn: (tt, 7r)-direction, S y = S x . Top row: Sequence of en- 
ergies io/uj res = 1.0 ... 1.4 at and above the resonance energy 
i^Jres — 0.5 J . Middle row: uj/cj res — 1.0... 0.7 at and be- 
low u> res . Bottom row: Energies far below w res , showing a 
crossover from parallel to diagonal incommensurability. Note 
the different vertical scale. 



A. Structure of the magnetic response in 
•wave-vector space 

Fig. M presents wave- vector scans of the magnetic sus- 
ceptibility Imx(q, w) for YBCO in the superconducting 
state at T = . For comparison with INS data an ex- 
perimental energy resolution 4r rs 5 meV has been sim- 
ulated through a quasi-particle damping T — 0.01J. 
For lu = uj res scans in parallel (tt, 0) as well as diag- 
onal (tt, 7r)-direction show a commensurate sharp peak. 
When energy is increased (top panel) , the resonance first 
evolves into a broad incommensurate structure with max- 
ima dispersing like spin waves P,36[ . At higher ener- 
gies around 2Ao « l.5u> res this turns into some feature- 
less background. When the energy is reduced from u> res 



10 



(middle panel), the peak also splits in parallel as well 
as diagonal direction, suggesting a circular structure in 
q-space. However, the maximum intensity is higher in 
(-7T, 0)-direction, reproducing the experimental observa- 
tion. The case u>/u> res « 0.7 can be compared to a study 
p5[ on underdoped YBa2Cu306.6 : The ratio of max. 
intensities I pa r/Idiag ~ 1-7 from Fig. H is consistent to 
the < 2.0 we read off the INS data in Ref. §| . The 
range 8 — 8 X — 0.052 . . . 0.065 of the displacement of 
peaks in {tt, 0)-direction is comparable to experimental 
values p3| reported [p5|,f37f for the same range of ener- 
gies uj/uj res = 0.8 . . . 0.7. In the bottom panel, where 
the energy is reduced even further, I par starts to weaken 
relative to Idiag ■ For uj/u> res < 0.35 the peaks in the 
diagonal {tt, 7r)-direction eventually dominate, as is ex- 
pected from a d-wave superconductor in mean-held the- 
ory at lu — > . However, such a crossover from parallel 
to diagonal incommensurability at low u> is not observed 
in experiment. At low energies the INS data indicate 
a strong isotropic suppression |9|,[37| , similar to what is 
seen |||| in LSCO . 

In the normal state, INS on the optimally doped com- 
pound shows a broad commensurate peak, its width de- 
pending weakly on energy pq . From the calculation at 
T > T c ss Td we get indeed a response almost indepen- 
dent of wave vector and energy, as is expected from the 
absence of nesting properties of the underlying Fermi sur- 
face (see Fig. @) . Note that both the commensurate (at 
io res ) and incommensurate intensity vanish in the normal 
state, as is observed fl37j in INS . The situation is differ- 
ent in underdoped YBCO , where the spin-gap regime 
is entered as T > T c . In mean-field theory the spin-gap 
phase is similar to the SC, with the SC gap becoming the 
spin gap (see Section III). Therefore the pattern in the 
magnetic response at u> < uj res persists at T > T c , which 
has also been observed in underdoped systems |3(],|80| . 
Experimental line shapes at T > T c are not reproduced 
in mean-field; however, we expect significant damping in 
the spin-gap phase if fluctuations are included. 



louin zone. Each d-wave node is surrounded by such a 
'banana-shaped' contour, and in an energy range around 
LUi nc rs 0.7uj res , where the parallel incommensurate pat- 
tern in Imx(q, Ui nc ) is most pronounced, these E(k) show 
almost flat pieces parallel to the magnetic zone boundary. 
This gives rise to a dynamic nesting contribution |8|,|8(| , 
which favors peaks at q = (tt ± 5, tt), {tt, tt ± 8) . In par- 
ticular, at lu = come a ratio of intensities I par / Idiag ^ 2 
is expected from the nesting argument, which is close to 
the value drawn from the numerical calculation (Fig. H) 
as well as experiment. An illustration of this effect has 
been given in Ref. Jl2| in Fig. 5 . 

In the preceding section it became apparent that the 
commensurate resonance at {tt, tt) depends on a suffi- 
ciently large n.n.n. hopping t' < . With respect to the 
(parallel) incommensurate pattern this is not the case: 
The underlying dynamic nesting effect is a general fea- 
ture of the d-wave superconductor. This is confirmed in 
a calculation of q-scans for t' = : The parallel incom- 
mensurability dominates for energies above the crossover 
from the diagonal one and below « 2\yJ | , where Imx 
becomes broad and commensurate. 
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FIG. 10. Wave-vector scans in {tt, 0)-direetion as in 

Fig. |, left column, but with 4r = 0.002 J . Shown are 
ui/uires = 1.0,0.9,0.5. Left: Bubble spectrum Imx" -7 " . 
Right: Magn. response Imx • The curves for 1.0 and 0.9 
are scaled xO.Ol ; the one for 1.0 features a 5-type peak, see 
text. 



B. The dynamic nesting effect 



An explanation of the incommensurate pattern below 
Lo res hi the superconducting state can be found in the 
dispersion -E(k) of the Gutzwiller-renormalized fermions 
|}42f . At vanishing energy to — > only particle-hole (ph) 
excitations 17 (q, k) with q connecting two d-wave nodes 
.E(k) > , E{k + q) > are possible. At finite hole dop- 
ing the nodes are shifted from (±7r/2, ±7t/2) towards the 
T-point (0, 0) ; thus the bubble spectrum Imx Jrr features 
peaks |82] at q = {tt ± 8', tt ± 8') , diagonally displaced 
from {tt, tt) . The curves for low co/uj res = 0.2, 0.3 in Fig. 
are still dominated by this type of ph-excitation. With 
increasing energy a different process gains importance, 
where q connects two contours E(k) — ui/2 in the Bril- 



C. Dispersion of the resonance 

Returning to the case t' = — 0.45i, for energies in- 
creasing towards the resonance energy uj res the incom- 
mensurate pattern eventually merges into the commen- 
surate resonance, as is seen in Fig. S, middle row. This 
is due to the final-state interaction Eq. (E4J) , which de- 
velops a pole at the commensurate position {tt, tt) for 
lj = cures ■ In Fig. no scans similar to Fig. O are made in 
{tt, 0)-direction, with the 'experimental damping' omit- 
ted. For uj/io res = 0.5 the bubble spectrum Imx"" r (q, w) 
features two humps in {tt, 0)-direction, and is zero outside 
these regions (left panel). This is due to a q-dependent 
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threshold $lo(q) = min/.0(q, k). With the interaction 
a J switched on, x trr ~ > X ? the humps are merely ampli- 
fied in intensity (right panel), and the nesting argument 
can be applied as above. When u> is increased, the two 
humps move closer. Additionally, sharp peaks appear in 
ImX > which eventually merge in a single 5-like resonance 
at (tt, 7r) for w —* co res . The incommensurate structure 
still present in the bubble spectrum for Lo/u res — 1.0 is 
completely superseded by the pole of Eq. ( pi[ ) . The lat- 
ter is driven by the bubble's real part Rex" r , which at 
U) = LOres is strongly peaked at (tt, tt) in q-space, assisted 
by the interaction J(q) = J[cos(q x )+cos(q y )] , which also 
favors the wave vector (tt, tt) . 
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FIG. 11. Dispersion in (tt, 0)-direction for x = 0.12. Pa- 
rameters as in Fig. H. Shown is a sequence of wave vectors 
&x — (0 . . . 0.12) 2-7T , S y = . Bottom panel: Bubble spec- 
trum Imx' rr ; Inset: Zoom of the threshold region. Middle: 
Inverse Stoner factor K x (— 1) . Top: Resulting response 
Imx from RPA . 

It is instructive to look at the commensurate- 
incommensurate crossover also in energy space, using 
the w-scans in Fig. O for several fixed wave vectors 
5 X = q x — tt , q y = tt . At q = (tt, tt) the bub- 
ble spectrum lmx irr (q, ui) features a single step-like on- 
set of spectral weight at the particle-hole (ph) excita- 
tion threshold $lo ; its consequences for the formation 
of the neutron resonance have been discussed in Sec- 
tion M above. In Fig. O (bottom) it is demonstrated 
how this threshold splits into two structures as we move 



away from (tt, tt) : the ph-threshold $7o(q) itself, shifting 
to lower energies, and a second step-like onset of addi- 
tional spectral weight at some fi 2 (q) > $7o(q) > which 
shifts up with 5 X . The denominator (real part) of Eq. 
© , K(q, W ) = [1 + o2 J[cos(q x ) + cos(q y )}Re X irr (q, lo)} 
thus shows a splitting of the corresponding log singu- 



larity into two (middle panel in Fig. 11), which in turn 
produces two peaks in the magnetic response Imx(q, w) 
(top panel). One of these peaks disperses to lower ener- 
gies and is identified as the q-dependent resonance, since 
it is located below the threshold $7o(q) an d is therefore 
sharp. With increasing 5 X its spectral weight is contin- 
uously reduced, since the height of the step at fio(q) m 
Imx IT,r decreases. The other peak near $72 (q) disperses 
to higher energies. It appears merely as a broad peak 
or shoulder within the damping continuum. It is respon- 
sible for the spin-wave like dispersion of broad maxima 
above u> res in Fig. |9| (top). The resonance, on the other 
hand, follows an upside-down dispersion and produces 
the incommensurate peaks below Lu res in Fig. |0j (middle 
row). 




FIG. 12. Dispersion of features of the bubble spec- 

trum Imx' rr and the spin-1 bound state (resonance) in Imx • 
8x = (Qx—tt) , q y — it . Shown is the ph-threshold Qo , the 2nd 
onset of weight f^2 , and the three v.H.s. associated with 2Ao . 
The position uj re .s of the resonance is indicated by squares. 
The ph-threshold f2o is also given in (tt, 7r)-direction as Qdiag ■ 

Results for the (tt, 0)-direction and doping level x = 
0.12 are summarized in Fig. O (top). It displays the 
ph-threshold Hq (q) , the second onset of spectral weight 
$^(q) , and the resonance position uj res (q) where it ex- 
ists. For comparison the ph-threshold fldiag(q) in (tt, t)- 
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direction is also shown. In contrary to firj(q) it has ze- 
roes at S x = ±{\f2kp — tt) , where q connects two d-wave 
nodes. In addition the energies of the van Hove singular- 
ities (v.H.s.) associated with 2An are given in the figure. 
At q = (tt, tt) the bubble spectrum Imx Jrr (Fig. 0) shows 
a single v.H.s. at W = 2A° w 0.72 J. From Fig. [ft] (bot- 
tom) it can be seen that this v.H.s. splits into three peaks 
with quite flat dispersion. 

The effect of strong underdoping is demonstrated in 
the bottom panel of Fig. n3 for i = 0.06 : The ph- 
threshold Slo at q = (tt, tt) and with it the resonance 
position shift down (compare also Fig. m. The q-range, 
where a sharp resonance exists, shrinks |57) , and the 
upside-down dispersion narrows. On the other hand, the 
v.H.s. around 2Ao as well as the maximum of 0,2 vary 
only weakly with doping (An increases slightly with un- 
derdoping). This will become important in the calcula- 
tion of wave-vector integrated susceptibilities in Section 
Vil below. 



D. Discussion 



The parallel type of incommensurability, i.e., a maxi- 
mum intensity at the points q = (tt, 7r±(J), (7r±<5, it) in the 
Brillouin zone is a generic feature of the d-wave SC state. 
In an energy range below 2A° the intensity is enhanced 
at these points due to the dynamic nesting mechanism 
|}42f . At very low energies, on the other hand, excita- 
tions from node to node in (it, tt) direction dominate |82| | 
and lead to a crossover to the diagonal type as to — » . 
It should be noted that the parallel incommensurability 
is not related to "stripes" |88|-|90[ , i.e., we do not con- 
sider the possibility of a combined ordering of charge and 
spin into quasi one-dimensional structures. Incommensu- 
rate pattern in Imx(q, iS) and their dispersion have been 
obtained with similar slave-particle methods |15 9l],|92| , 
BCS theory ||||,||L or the FLEX approximation for 
the Hubbard model [20,p4| . The present slave-particle 
approach predicts how the dispersion of the resonance 
(associated with f2o(q) m Fig. [12h and the spin- wave like 
dispersion of broad peaks above it (following ss 1^2 (q)) 
change with underdoping. When x is reduced, the peaks 
connected to ^(q) should be observable in a wider en- 
ergy range above the resonance. Furthermore, near the 
bottom of f^2 at q = (tt, tt) the density of states and thus 
the damping is reduced, leading to sharper peaks in the 
underdoped case. Experiments [§J3J] actually indicate 
that dispersing 'spin-wave' peaks above uj res are better 
resolved in the more underdoped sample. 



VII. EFFECT OF THE COUPLING IN THE 
DOUBLE-LAYER 

So far we have considered a single Cu02 layer as the 
most important structural element of cuprate supercon- 
ductors. However, YBCO and BSCCO contain two cou- 
pled CuC>2 planes per unit cell. The observed suscepti- 
bility actually follows §HI3@ 

X(c[,q z ,uj) = x + (q,^)cos 2 (-g z ) +x _ (q,^)sin 2 (- gz ) 

This form reminds of the odd (— ) and even (+) linear 
combination of spin waves in the undoped parent com- 
pound (9^^3l ■ 1 ^ s t ne in-plane wave vector as before, d 
is the distance of the planes within a bi-layer sandwich. 
We used the single-layer model as an effective model for 
the odd ("acoustic") susceptibility \~(q,oj) . The exper- 
imentally observed neutron spectra in the odd mode, in 
particular the resonance at q = (tt, tt) and its doping de- 
pendence are well reproduced by the single-layer model. 
Our description of the resonance in Inr^ - ((tt, tt), to) does 
not rely on the bi-layer structure of the material. The 
important ingredient is the topology of the underlying 
Fermi surface in combination with the d-wave super- 
conducting state. This has been discussed in detail in 
Section M above. The even ("optical") mode spectrum 
1vox + ((tt, tt), lu) appears different in experiment. It shows 
merely a broad peak with dim intensity [pf . In this sec- 
tion the calculation is extended to the bi-layer system. It 
is shown that the suppression of the resonance in the even 
mode is mainly a consequence of the inter-plane exchange 
coupling J^- . The odd-mode susceptibility, on the other 
hand, resembles the one obtained from the single-layer 
model. 

The bi-layer modes have also been explored by averag- 
ing the experimentally measured magnetic response over 
the in-plane Brillouin zone fip4,|dLp5| , 



Im %2DH 



(2TT) 



Imx (q,w) 



(29) 



After the q-integration has been performed the odd mode 
spectrum Imx^^u;) is still dominated by a sharp reso- 
nance; it occurs at the same energy uj res < 40 meV as in 
Imx~ ((tt,tt),u)) , but with significantly diminished am- 
plitude H . In the even (+) mode a second energy scale 
becomes apparent. Iul^d shows no resonance, but a 
rather broad peak pUo or soft onset of spectral weight 
p4j . The location uihump ~ 80meV of this 'hump'-like 
structure is almost independent of doping, in contrast to 
the strongly doping dependent U) rea . The 'hump' also 
appears more El[ or less M clearly in the odd (— ) mode. 
Xid wm b e studied later in this section. It turns out that 
particle-hole (ph) excitations across the maximum gap 
A lead to a hump-like peak in both modes in Imx 2D , 
at an energy % 2A° almost independent of doping. 
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A. Results for the bi- layer system 

Theoretical expressions for the susceptibility of two 
coupled planes have been derived in Section n. From 
Eqs. (p2,pQ) the mode susceptibilities are given by 



x ± (q,^) = 



xi rr M 



i + j ± (q)xir r M 



(30) 



in units of (fl^s) 2 • xZ r is obtained from Eq. (|23| ) with 
p = (q, p z ) and p z = {0, n} for the modes {+, — } . The 

X differ in their respective irreducible particle-hole bub- 
ble x l p T an d the effective interaction Eq. (g6|) , 

J^q) = aJ(q) ± J X , J(q) = 2J[cos(q x ) + cos{q y )} 

For the in-plane parameters we take a = 0.35, t — 2 J, t' — 
— 0.45t as before, and for the coupling of the two CuC>2- 
planes within a bi-layer we chose an antifcrromagnetic 
exchange J 1 - = 0.2 J and an inter-plane hopping [{45| , fl6| 

t-'-(q) = 2t ± [cos(q x ) - cos(q y )} 2 + t& 

with t x = O.lt and ^ = 0. 

We assume an in-plane superconducting order param- 
eter A with equal amplitude and phase in both layers. 
The self-consistent solution of the mean-field equations 
( |l9| ) then leads to a vanishing inter-plane gap A ±0 = 
(Ref. Eqs. (|||Jl|)) and a very small uRVB amplitude 

X I ~ 0.03 , which has been defined in Eq. @ . There- 
fore the influence of t , J on the fermions that consti- 
tute x l p T is merely a small splitting of the bandstructure 
into bonding and anti-bonding bands through the effec- 
tive inter-plane hopping t 1 - (k) m xt (k) . 

40 

c^ 3 30 

CD 

ffl 20 

ZS 
O.10 

E 
< 


0.2 0.4 0.6 0.2 0.4 0.6 

Energy / J Energy / J 

FIG. 13. Odd- and even-mode susceptibility (imag. part) 
of a bi-layer system, for fixed in-plane wave vector q = (n, n) . 
Parameters are t = 2J, t' = -0.45t, a = 0.35, T = 0, 
4r = 0.04 J « 5 meV as in the preceding sections, and 
J± =0.2J, tx = 0.1*. 
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Results for fixed in-plane wave vector q 



(7T,7r) 



are presented in Fig. |13l , with 'experimental' energy 
resolution 4r = 5 meV . A resonance appears in the 
odd mode susceptibility, which varies with doping as 
in the single-layer case. The even mode, on the other 



hand, shows a broad peak with much reduced intensity. 
This is mainly due to the mode-dependent interaction 
|J + (7r,7r)| < |J-(7r,7r)| in Eq. (§(]), which shifts the 
pole in the even (+) mode into the damping continuum. 
The damping effect is supported by the above-mentioned 
splitting of fcrmion bands. It should be emphasized that 
the resonance in the even mode is not totally suppressed, 
but shifted and strongly damped. In experiment B a 
similar observation has been made. 
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FIG. 14. Wave-vector q integrated odd- and even-mode 
susceptibilities IniX2D from Eq. (J3lJ) for hole filling 

Inset: 

The max- 



from Eq. @ for 
x — 0.06... 0.18. Parameters as in Fig. 
q-integrated bubble spectrum Imx^D f° r X = 0. 



imum is located close to 2A 



o 



0.78J. 



A new feature appears if we look at the wave-vector 
integrated susceptibility Eq. (E9) , which is shown in Fig. 
H . Similar to the case of fixed q = (tt,?!") a resonance 
appears only in the odd mode. It appears at the same 
position as in Imx~ (tt, tt) with the same strong doping 
dependence. Additionally, both modes Iniy^u^) show 
a broad peak ('hump') at an energy somewhat below 
2A° , almost independent of doping (2A° sw 0.78 J for 
x = 0.08) . For an explanation of this 'hump' we first go 
back to the single-layer case: The spectrum lmx"" r (q, oj) 
of the irreducible ph-bubble shown in Fig. O (bottom) is 
dominated by peaks around 2A° . These van Hove singu- 
larities (v.H.s.) follow the quite flat q-space dispersion 
shown in Fig. O. When the wave- vector is integrated 
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2 

over in Imx^oCw) = JJ w^p I m x' rr (qi w ) > the v.H.s. con- 

tribute a large density of states, leading to a broad peak 
with maximum at w = 2A° . This argument extends to 
the bi-layer system: The hump appears almost identi- 
cally in both modes of the bubble Imx^D^i 1 - ) i which is 



shown in the inset of Fig. 14 . Its position follows 2A° to 



slightly higher energies for reduced doping level x . 

When the final-state interaction Eq. (^) is switched 
on in the odd (— ) mode, the resonance appears (see Fig. 
[14 top). Since spectral weight is shifted to lower ener- 
gies, the hump is relocated to an x-independent position 
below 2A° . The even (+) mode (Fig. pTl bottom) experi- 
ences a weaker renormalization through Eq. (pG) , no res- 
onance is formed, and its hump is relocated less strongly. 
Recent FLEX calculations for a 3-band single-layer Hub- 
bard model in the overdoped regime p6| give results for 
ImX2£>(w) comparable to our odd mode susceptibility. 

The intensity of the resonance in Imx^ D ( w ) IS much 
reduced from its value in Imx~((7r, tt),uj) measured at 
fixed wave vector. Therefore in Im^D the resonance 
at uj res and the excitations across the maximum gap 
at tOhump ^ 2A° are of comparable intensity and can 
both be observed experimentally. This is due to the fact 
that the latter occupy a large part of the in-plane Bril- 
louin zone (BZ), while the resonance is just a narrow 
peak in q-space. Despite its large amplitude, the reso- 
nance contributes only little to Eq. (J2^) . For the actual 
computation of Imx 2 D we use d wave-vector scans |9S|] 
in (n, 0)-direction like those shown in Fig. || (extended to 
the whole BZ) and assumed [Q a susceptibility isotropic 
around (n, n) , i.e., 



X^b( w ) 



2tt 



fcdfcx ± (q, w) , q = (t + k,ir) (31) 



The resonance is actually so sharp in q-space (see 
Fig. fj) that it does not become visible in the result- 
ing ImX2x>(w) ■ For each u> the respective q-scan has 
therefore been convoluted with a Gaussian of FWHM= 
0.25r.l.u.= 0.57T, in order to simulate a finite instrumen- 
tal wave-vector resolution. Application of Eq. (plh then 
leads to the curves shown in Fig. pi . 



B. Comparison to experiment 

Two experimental groups studied the wave- vector inte- 
grated magnetic response Imx 2 D in underdoped YBCO . 
Refs. Hj9q] reported a lineshape for YBCOg.6 which 
agrees quite well with the theoretical result Fig. H4 
for x < 0.08 . A 'hump' in ImxJo (even) appears at 
R3 100 meV, I111X2D (odd) shows a similar structure at 
a somewhat lower energy w 90 meV . The well-known 
resonance appears only in Im^,^ , at 34meV. In Refs. 
P,p4l two underdoped samples YBC0 6 . 7 and YBC0 6 . 5 



have been studied. In the even ("optical") mode of 
YBCO6.7 a hump appears around 70meV, whereas the 
odd ( "acoustic" ) mode shows a weak hump-like structure 
at rs 55 meV , separated from the resonance at 33 meV . 
In the more underdoped sample YBCO6.5 these features 
tend to move up in energy, while the resonance in Imx^ 
shifts down to 25 meV . 

Although the detailed experimental lineshapes are not 



unique [ 100 1 , the qualitative features of our calculation 



are found in the neutron-scattering spectra. In particular 
we obtain the different dependence on doping level of the 
resonance at u> res in the odd mode and the hump-like 



hump 



in both modes. 



A1S0 1S W Wp ° f the 



odd mode lower than the uj^ u of the even mode. The- 
ory and experiments can also partly be compared quan- 
titatively. The measured neutron-scattering intensities 
P) Sj| a re of the same order as the theoretical ones in 
Fig. |l4j (using J = 120meV, i.e., l/i|/J = 8.3/i|/eV). 
The maximum of the hump in the even, odd mode in Fig. 
HI occurs at u> + ~ w 0.6 J, 0.53 J = 72meV, 64meV, in 
good agreement with the measurements || on YBCO6.7 
at low temperature. Note that the maximum gap A = 
30 — 45 meV M is consistent to the value from the mean- 
field calculation (see Section |H|) , 2A° is the upper limit 
for the hump position in Fig. |14|. 



C. Discussion 

The wave-vector integrated magnetic response in un- 
derdoped systems is characterized by two energy scales 
with opposite dependence on the doping level. The first 
is the position uj res of the "41 meV resonance" , which 
appears in the odd (— ) mode for fixed wave vector 
q = (-7T, it) as well as in the q-integrated susceptibility. It 
moves down in energy when doping x is reduced and be- 
comes a Bragg peak at the transition to the AF ordered 
state at x — x c . The second is essentially the maxi- 
mum gap A , which increases with reduced x . It deter- 
mines the position UJhump °f the additional broad peak 
('hump'). The latter appears in both bi-layer modes, but 
only if the in-plane wave- vector is integrated over. The 
hump is caused by particle-hole excitations across 2A° , 
and is pulled down somewhat by the final-state interac- 
tion Eq. (|3(j) . It should be noted that it is very robust 
against a variation of the next-nearest-neighbor hopping 
t' , i.e., the topology of the fermion's band structure and 
Fermi surface. Whereas the resonance vanishes for t' = 
the hump remains almost unaffected. The mechanism 
is very much different from the optical spin waves that 
appear in the undoped x — bi-layer system if a finite 
Neel order parameter is taken into account. Therefore 
the appealing similarity of the hump- (or threshold-like) 
feature in the superconducting YBCO samples and the 
optical spin-wave gap seen in the undoped parent com- 
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pound |97j is accidental. 



VIII. SUMMARY AND OUTLOOK 

This paper presented a theory for the magnetic 
excitation spectrum of YBa2Cu306+i, (YBCO) and 
Bi 2 Sr 2 CaCu208+,5 (BSCCO) superconductors. We con- 
sidered the so-called "41 meV resonance" at fixed in- 
plane wave vector q = (ir, it) , the magnetic response in 
q-space, the pecularities due to the bi-layer structure of 
YBCO and BSCCO, and the local (q-integrated) sus- 
ceptibility. Most of the results are in good agreement 
with the neutron-scattering experiments. The resonance 
is obtained as a collective spin-1 excitation in the super- 
conducting and spin-gap states (the latter corresponding 
to the pseudo-gap regime of cuprates). Its energy scale 
and spectral weight as function of the doping level x at 
low temperature are satisfactorily reproduced. The ab- 
sence of damping in optimally doped systems is caused 
by the d-wave superconducting gap in connection with 
the hole-type topology of the underlying Fermi surface. 
The bi-layer structure is not necessary for the resonance 
to form in the odd-mode (acoustic) susceptibility. The 
mere effect of the finite inter-layer coupling J 1 - is an al- 
most suppression of the resonance in the even (optical) 
mode. The observed pattern of incommensurate peaks in 
q-space has been traced back to a dynamic nesting effect 
of the d-wave superconductor, and the peak's dispersion 
has been derived for optimally and underdoped systems. 
Besides the resonance a second, hump-like feature ap- 
pears in the wave-vector integrated magnetic spectrum. 
It is caused by particle-hole excitations across the maxi- 
mum gap A that occupy a large area in the 2D Brillouin 
zone. Their energy oJhump ^ 2A° is almost independent 
of hole filling, in strong contrast to the resonance position 
to res . Also is this hump insensitive to the Fermi-surface 
topology. 

A salient property of the resonance at (tt, tt) is the 
variation of its energy uj res with hole filling x . The t—J- 
model, i.e., the doped Mott insulator naturally provides 
the energy scales for the resonance in the underdoped 
and overdoped regimes: The mean-field theory describes 
magnetic excitations in terms of quasi particles (QP) (the 
fermions) with a Fermi velocity Vf ~ {x + 0.15J/t)vF 
reduced from the bare parameter. Hence in underdoped 
systems the QP's chemical potential is smaller than the 
gap, \n^\ < A , and determines the scale for the reso- 
nance energy as to res < 2\/j,f\ (see Fig. [7p. Thus uj res is 
found to increase with hole filling x , in accordance with 
experiment. In the overdoped regime, on the other hand, 
we have \n*\ » A , and the resonance uj res < 2A° is 
connected to the gap which decreases with x . 

The mean-field theory in its present form overestimates 
the antiferromagnetic (AF) state in the phase diagram. 



Therefore we had to introduce the phenomenological pa- 
rameter a , which reduces the interaction J — ► a J in 
the spin-flip particle-hole channel Eq. (|2J) of the quasi 
particles. The present study shows that already the sim- 
plest model a(q, u), x) — a leads to consistent results 
for neutron-scattering spectra and magnetic correlation 
length £af (x) in the relevant range of doping, energy, and 
wave vector. With respect to the doping dependence of 
Lu res and £af this is due to the above-mentioned renor- 
malization of the QP . In the half-filled case x = , which 
we did not consider here, the mean-field theory delivers 
a Neel state with the correct spin- wave velocity only if J 
is kept unrenormalized, i.e., a(x — 0) = 1 . With doping 
the AF state is destroyed by the propagation of holes in 
the spin background fl6q,[71 72 . We expect that a refined 
theory, where these processes are included as corrections 
to mean field, yields an a(x) which descreases quickly in 
the AF region x < x c and then levels off in the param- 
agnetic (superconducting) phase x > x c . This is subject 
to future work. 
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